Large sample-to-sample fluctuations of the nonequilibrium critical current through 

mesoscopic Josephson junctions 
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We present a theory for the nonequihbrium current in a mesoscopic Josephson junction which 
is coupled to a normal electron reservoir, and apply it to a chaotic junction. Large sample-to- 
sample fluctuations of the critical current 7c are found, with rms Ic — y/NeA/h, when the voltage 
difference eV between the electron reservoir and the junction exceeds the superconducting gap A 
and the number of modes A'^ connecting the junction to the superconducting electrodes is large. 

PACS: 74.50.+r, 74.20.Fg, 74.80.Fp 

Over the last years there has been an increased 
interest in the nonequilibrium Josephson current in 
mesoscopic multiterminal superconductor-normal metal- 
superconductor (SNS) junctions. Nonequilibrium in the 
junction is created by quasiparticle injection from one or 
several normal electron reservoirs, connected to the nor- 
mal part of the SNS junction. By controlling the voltage 
applied between the normal reservoirs and the SNS junc- 
tion, it has been shown in recent experiments that the 
Josephson current can be suppressed [Q,^, reversed |^] 
and in the case with injection from a superconducting 
reservoir, even enhanced jj]. 

The microscopic mechanism for these effects, nonequi- 
librium population of the current-carrying Andreev lev- 
els, was discussed by van Wees et al. Q already in 1991. 
Thereafter, the nonequilibrium Josephson current in var- 
ious multiterminal geometries has been studied in both 
diffusive |^-^ and quantum ballistic junctions. In 

Ref. pO[ | it was pointed out that the nonequilibrium 
Josephson current in ballistic SNS junctions can not be 
described only in terms of the nonequilibrium population 
of Andreev levels: There is also a quantum-interference 
addition to the Josephson current, which results from the 
difference between the scattering-state wavefunctions for 
injected electrons and holes. 

In this paper we develop a general theory of the 
nonequilibrium Josephson current in three-terminal SNS 
junctions (see Fig. [I]), in terms of the scattering matrix 
for electrons and holes injected from the normal reser- 
voir. The theory is then applied to a chaotic junction, in 
the limit of weak coupling to the normal reservoir and at 
zero temperature. We find that the quantum-interference 
contribution gives rise to sample-to-sample fluctuations 
of the critical current Ic which are much larger than the 
equilibrium fluctuations ||l^,|l^: For a large voltage V 
(with eV A, the superconducting gap), 

peA 



FIG. 1. Three-terminal SNS junction, consisting of a 
mesoscopic scatterer (grey shaded) connected to two super- 
conducting reservoirs via contacts 2 and 3 and a normal reser- 
voir via contact 1. The black bar in contact 1 indicates a 
tunnel barrier, the arrows the direction of positive current 
flow. 
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hence the fluctuations are of the order of the ensemble- 
averaged critical current itself. (Here N is the number of 



modes connecting the junction to each of the supercon- 
ducting electrodes.) In this regime the current results 
from the quantum-interference contribution alone. For 
eV < A the critical current is of order N^, with fluc- 
tuations of order 

n 

A model of the junction is presented in Fig. n^. A 
mesoscopic scatterer is connected to two superconducting 
leads via ballistic contacts, each supporting N transverse 
modes. The phase difference between the superconduc- 
tors is 0. The scatterer is also connected to a normal 
reservoir via a contact with M modes, containing a tun- 
nel barrier with transparency F. A voltage V is applied 
between the SNS junction and the normal reservoir. We 
assume that the resistance of the injection contact is the 
dominating resistance of the junction, such that the po- 
tential drops completely over the injection point. In or- 
der to preserve nonequilibrium, the strength of the tunnel 
barrier F is however limited by the requirement that the 
dwell time of the injected quasiparticles cx 1/F must be 
smaller than the inelastic scattering time in the junction. 

Under these conditions, the distribution of the quasi- 
particles in the junction is determined by the distribu- 
tions n°^'^-' = riF(£' T eF) of electrons (holes) in the 
reservoir, where n-p — [1 + exp{E/kT)]^^ . The current 
in contact j = 1, 2, 3 can then be written as 



1 



(2) 



with 



the current density of the scattering states re- 



sulting from injected electron (hole) quasiparticles from 
the normal reservoir and i^- the total current density for 
quasiparticles injected from the superconductors {V^ = 
for subgap energies \E\ < A). 

The current I = {I2 + -^3)/2 flowing between the su- 
perconductors can be rewritten by using the current con- 
servation for each energy i^'^ + 12'^ — ig and the fact 
that no current is flowing in the injection lead in equilib- 
rium, i<l+i\ = 0. It takes then the form / = 1"""^ + r'i, 
where the equilibrium current at = is given by 
r'i = J dE[i+ + [il + i|)/2]nF, and 



rnoq _ 



dE 



y(n° + n'^-2nF) + y(n^-nh) 



(3) 



Here the current densities 



= {^2 



i\ -{■ 12 



iq and 



13)72 are the sum and the dif- 



ference of the current densities of the scattering states 
for injected electrons and holes. The contribution cx z+ 
to results from the noncquilibrium population of 
the Andreev levels, while the current oc i~ accounts 
for the quantum-interference contribution as well as for 
an asymmetric splitting of the injected current Ii = 
] dE{i\ - i\)(n° - n^) /2. 

We will now express the current densities in terms 
of the scattering matrix S of injected quasiparticles 
from the reservoir pl] |. The current densities are cal- 
culated most conveniently in the contacts j — 1, 2, 3, 
where the wavefunctions are plane-wave solutions to the 
Bogoliubov-de Gennes equation. A wave incident on the 
scatterer from leads 2 and 3 is described by the AN vector 
of wavefunction coefficients cin — (cj ^, Cg'", Cj'", Cg'^). 
The superscript ) denotes a positive (negative) sign 
of the wave vector. Correspondingly the outgoing wave 
is given by Cout = (C2' , Cg'^, C2^ , ^')- At the NS inter- 
faces, Andreev reflection is described by the scattering 
matrix 



Sa = a 
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(4) 



such that Cin = •S'ACout, with a — exp[— i arccos(£'/A)]. 
The wavefunctions in the three contacts are then 
matched with help of the {2N + M) x (2N + M) scatter- 
ing matrix S' of the normal region (including the tunnel 
barrier). 
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^'M, (6) 
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such that Cout = S'NCin, and matrices which involve also 
contact 1, 



T = 



ti2{E) h3{E) 
tU-E) tU^E) 

(t2l{E) 





hi{E) 
, tl^(-E) 
V tl^(-E) 



n = 



rii{E) 
Tl^{~E) 



^From these ingredients, the coefficients c can be calcu- 
lated and the current densities in Eq. (^ are obtained 
from the quantum mechanical expression for current. 
The current densities i+ and i~ follow after some ma- 
trix algebra, 

i+{E) = ^tr (r't(l - 44)-V;(l - SnS^r'T') , 
i-{E) - ^tr (T't(l - 5l5i)-V;(l - 5N5A)-^T'a,) , 

(7) 

where CTz = diag (1, 1, —1, —1), a'^ — diag (—1, 1, 1, — 1). 
The scattering matrix S for injected quasiparticles from 
the normal reservoir can be written as 



S = 



^cc ^hc 



(8) 



After some further matrix algebra we obtain from Eqs. 
(0) and (^ the expressions (for subgap energies \E\ < A) 



in 



^tr (S^isr.^ 
in 



(9) 



with Tz = diag (1,-1). (The expression for i'^ is well 
known [||jl4|.) Eqs. (|) and (|) are our general results 
for the noncquilibrium Josephson current. 

In general, the current flowing between the supercon- 
ductors contains also the part of the injected current 
which is asymmetrically split between contacts 2 and 
3 This is not the case when the SNS junction is 

weakly coupled to the reservoir (F ^ 1), because the in- 
jected current is then negligible. In this limit the matrix 
5n = S^Q -\- F(5S'n can be expanded to first order in F, 
where iSno is unitary. The two current densities 1+ and 
i~ have the same discrete spectrum of Andreev levels, 
given by the solutions E^ of det(l — S'a-S'no) = 0, but dif- 
ferent spectral weights. The current density i"*" reduces 
to the well-known expression for the closed junction, 



-{E)=Y,I+5{E-E,,), 



4- = ¥^- (10) 



We introduce a non-unitary matrix S'n, describing only 
the scattering between the contacts j = 2 and 3, 



The current density i can be found from first-order per- 
turbation theory in the tunnel-barrier transparency F, 
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FIG. 2. Ensemble-averaged critical current (7c) (solid thick 
line) and the fluctuations rms 7c (dashed thick line) as a func- 
tion of voltage V between the normal reservoir and the junc- 
tion. The thin lines are the result with 7^ = in Eq. (p^. 
The junction has = 10 modes to each of the supercon- 
ducting electrodes and A7 — 1 mode to the normal reservoir. 
Inset: the voltage range 0.965A <eV < I.OIA. (10^ random 
matrices S' have been generated). 



(E) = ^ I-S{E - i?„), /- = i?„/+, (11a) 
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where the unitary matrix U diagonaUzes the unitary ma- 
trix product S'a'S'no = Udia,g{X)U^ . One can show with 
help of the corresponding eigenvalue equation that the 
ratios < 1. It should be pointed out that the matrix 
(55n can not be expressed in terms of the closed junction 
scattering matrix S'no, i- e. the current density i~ depends 
manifestly on the properties of the contact between the 
normal reservoir and the SNS junction. 

In order to investigate the mesoscopic fluctuations of 
the nonequilibrium current in more detail we now ap- 
ply our theory to a chaotic SNS junction^ in the limit 



of weak coupling to the normal reservoir 1 14 1 . Here we 
only consider the simplest case, in which the dwell time 
in the normal scatterer (with the superconducting leads 
replaced by normal ones) tdwcU < ?i/A. (Our main con- 
clusions should apply also for the opposite case.) For 
such a junction we can neglect the energy dependence of 
S", which is then distributed with the so-called Poisson 
kernel P(S") cx |det(l-(S"t)S")|"(2A'+*^+i), where {S') is 
the ensemble-averaged scattering matrix jlGf . (The mag- 
netic field B = 0, which gives a symmetric scattering 
matrix 5" = S''^.) Furthermore, the current density for 
energies outside the gap vanishes [|ll|. Using the energy 
symmetries i^{E) — —i^{—E) and i^{E) = i^{—E), the 
total current at zero temperature. 



/ = - 



E 

E^>eV 



E„<eV 



r 



(12) 



can be written as a sum over the currents and /~ 
carried by the individual Andreev levels with positive 
energies i?„. Eq. (|l2|) provides a simple picture where 



in equilibrium all Andreev levels carry the currents /+. 
Increasing the voltage, the Andreev levels one by one 
switch from /+ to /~ when the voltage is passing through 
eV = En- At eV > A, all levels carry the current I~ . 

In terms of the transmission eigenvalues T„ of the ma- 
trix 5*0, the Andreev bound-state energies are given by 
En = A(l - r„ sin^ 0/2)^/^, hence the relation 



/+ = -(eA/2;i)r„ sin0(l - r„ sin^ 0/2) 



-1/2 



(13) 



The statistical properties of the equilibrium current I'^'^ — 
-En^rt are known Q, with (J'^i) ~ 7VeA/?i and 
rms I^"^ ~ e/S./fi. 

For eV" > A the current is / = ^r7 = Z^n^n^n ■ 
The statistics of the ratios i?„ follows from the construc- 
tion of all perturbations 5S^ which are compatible with a 
given S'no (i-e., both matrices follow from the same scat- 
tering matrix of the open scatterer ) . For M = 1 such 
an analysis results in 

Rn = (1 -r„)i/2[(sin2 0/2)-i -r„]-i/2sin/3„, (14) 

where the angles (parameterizing the coupling to 

the reservoir) are independent random numbers, uni- 
formly distributed in the interval [0,27r). As a conse- 
quence for fixed phase difference (j) the average current 
(/) = 0, and the fluctuations rms / ~ ^/NeA/h because 
is a sum of A'^ independently fluctuating numbers 
The precise value of the fluctuations can be calcu- 
lated upon replacing the sum in {P) = {J^ni^riln)^) 
(valid due to the independence of the /3„) by an in- 
tegral over the transmission eigenvalues, with density 
p{T) = 7V7r-i[T(l - T)]-i/2. This results in 



rms / 



NeA 



!sin2 0/4 sin^(f)/2 
o 



2?i|tan0/2|y 2 COS0/2 COS0/2 
for eV > A, (15) 



which is parametrically larger than the equilibrium fluc- 
tuations when ^ 1. 

Another physical quantity of interest is the critical cur- 
rent Ic , the largest possible current for a given realization. 
Because of /(</)) = —I{—4') it sometimes makes sense to 
restrict the phase to < ^ < tt and to consider the cur- 
rent which is largest in modulus; Ic can then be positive 
or negative |^ . (With this definition, the average critical 
current vanishes for eV > A.) In the following, how- 
ever, we maximize over — tt < < tt, hence Ic is always 
positive, as it is obtained from the I /V characteristic in 
experiments. The ensemble-averaged critical current and 
its fluctuations (obtained from a numerical simulation of 
the random-matrix ensemble with = 10 and M ~ 1) 
are shown in Fig. H, as a function of applied voltage eV. 
The result is compared to the contribution of /"'" in Eq. 
(|l^ ) alone, which only takes the nonequilibrium popula- 
tion of the Andreev levels into account. 

For < eV ^ 0.54A the critical current is equal to its 
equilibrium value, because at the non-fluctuating criti- 
cal phase [O] — 2 all bound-state energies En > eV 
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FIG. 3. Ensemble-averaged critical current (7c) (open cir- 
cles) and the fluctuations rms 7c (full circles) as a function of 
the number of junction modes A'' for a single injection mode 
M = 1 (upper panel) and the number of injection modes M 
for A'^ = 10 junction modes (lower panel). The curves are 
oc N^^'^ (upper panel) and oc M~^l'^ (lower panel). 



(in general the energies lie in the interval [Acos(^/2, A]). 
In the range 0.54A < eV < 0.98A the critical phase 
is determined by the condition cos^c/S = eV j IS. that 
the first Andreev bound state drops below el^, with only 
small fluctuations due to the high density of transmis- 
sion eigenvalues r„ « 1. Hence the critical current is 
Jc = I'^'^{(l>c)- In this regime the quantum- interference 
contribution /~ in Eq. ( |l2| ) does not play any role be- 
cause (I"^) ^ rms I^. For a voltage eV « 0.98A very 
close to the gap, and /~ are both of order \/]VeA/?i, 
and the critical current starts to deviate from what one 
would expect from a pure nonequilibrium population of 
the Andreev levels. (For increasing N the cross-over volt- 
age eV — > A.) In parallel the fluctuations of the critical 
current increase. The critical current remains constant 
for eV > A, where it is given solely by /~. 

The critical current for eV > A and its fluctuations 
as a function of junction modes N are shown in the 
upper panel of Fig. ^ The mean critical current is 
(/c) ~ 0.16y/NeA/h. The fluctuations arc of the same 
order, rms Ic — 0.1\lNeA/h, which is by a factor of 
about ^/N /3 larger than the equilibrium fluctuations. 
Hence the N dependence in Eq. ( |l5| ) carries over to the 
average critical current and its fluctuations. 

Finally let us consider the dependence of the critical 
current on the number of injection modes M . This num- 



ber is significant because the current /~ depends mani- 
festly on the coupling of the reservoir to the junction [see 
Eq. (0)], in contrast to the current which only de- 
pends on properties of the decoupled junction. The lower 
panel of Fig. || shows that the critical current and its fluc- 
tuations at eV > A are suppressed when M is increased. 
The functional dependence is approximately cx AI^-^^^. 
The curves flatten out when M becomes larger than the 
total number 2N of modes connected to the supercon- 
ductors. Thus, for an experimental observation of the 
large fluctuations predicted above, an injection contact 
with few modes is favorable. 

In conclusion, we have studied the nonequilibrium 
Josephson current in a mesoscopic SNS junction con- 
nected to a normal electron reservoir. It is found that 
the current can be expressed in terms of the scattering 
matrix for the quasiparticles injected from the normal 
reservoir, Eqs. (^ and (^. As an application we con- 
sidered the nonequilibrium current in a chaotic Joseph- 
son junction at zero temperature, weakly coupled to the 
normal reservoir. It is found that the fluctuations of 
the critical current for a voltage eV > A are of order 
rms Jc ~ ^/NAe/h, which is of the same order as the 
mean critical current itself, and much larger than the 
equilibrium fluctuations (of order Ae/h). 
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